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The standard generic quantum computer model is studied analytically and numerically and the 
border for emergence of quantum chaos, induced by imperfections and residual inter-qubit couplings, 
is determined. This phenomenon appears in an isolated quantum computer without any external 
decoherence. The onset of quantum chaos leads to quantum computer hardware melting, strong 
quantum entropy growth and destruction of computer operability. The time scales for development 
of quantum chaos and ergodicity are determined. In spite the fact that this phenomenon is rather 
dangerous for quantum computing it is shown that the quantum chaos border for inter-qubit coupling 
is exponentially larger than the energy level spacing between quantum computer eigenstates and 
drops only linearly with the number of qubits n. As a result the ideal multi-qubit structure of the 
computer remains rather robust against imperfections. This opens a broad parameter region for a 
possible realization of quantum computer. The obtained results are related to the recent studies of 
quantum chaos in such many-body systems as nuclei, complex atoms and molecules, finite Fermi 
systems and quantum spin glass shards which are also reviewed in the paper. 
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I. INTRODUCTION 

On the border between two Millennia it is natural to 
ask a question, what will be the origin of future human 
power? Even thirty or twenty years ago the standard an- 
swer would be: nuclear. But now in a view of amazing 
computer development all over the world it becomes clear 
that the future power will be related to ability to count as 
fast as possible. In I Millennium this ability was basically 
comparable with finger counting, while at the end of II 
Millennium it made enormous jump with computer cre- 
ation which led to a qualitative change in human society. 
During last two decades the power of modern computers 
demonstrated a constant impressive growth due to tech- 
nological progress and creation of chips of smaller and 
smaller size. In a near future this size should reach a 
scale at which the quantum nature of physical laws will 
become dominant. As a result, we unavoidably come to 
the creation problem of quantum computer. Such a com- 
puter should be essentially based on quantum mechanics 
and operate with unitary transformations and quantum 
logic. The unitary nature of transformations allows to 
exclude energy dissipation that should play an impor- 
tant role on small scales. At the same time, as stressed 
by Feynman p| , the classical computer has enormous dif- 
ficulty in simulation of many-body quantum systems due 
to exponential growth of the Hilbert space with the num- 
ber of particles and hence, of the computational efforts. 
Due to that it is possible to expect that a computer com- 
posed of quantum elements will be much more efficient 
for solution of quantum, and may be other, problems. At 
present quantum computer is viewed as a system of cou- 
pled n qubits being two-level quantum systems or one- 



half spins (see recent review and references there in). 
The computer operation is based on controlled series of 
two-qubit coupling switch on and off which together with 
one-qubit rotations allow to realize any unitary operation 
in the Hilbert space of size Njj = 2 n . In this respect 
the inter-qubit coupling becomes unavoidable property of 
quantum computer. 

Recently a great increase of interest to quantum com- 
puting has been generated by the work of Shor p] who 
constructed a quantum algorithm which performs large 
number factorization into primes exponentially faster 
than any known classical algorithm. Also Grover showed 
|^| that a search of an item in a long list is done much 
faster by quantum computer. The enormous gain in 
the computation rate is reached due to high parallelism 
of multi-qubit quantum evolution and quantum interfer- 
ence. Together with a recent theoretical development of 
quantum error-correcting codes these exciting re- 
sults stimulated various experimental proposals for re- 
alization of quantum computer. The variety of phys- 
ical systems proposed is really amazing and includes: 
ion traps Ml, nuclear magnetic resonance systems [ [To| , 
nuclear spins with interaction controlled electronically 
|pd| , ^2| or by laser pulses Jl3| , quantum dots JlJ] , Cooper 
pair boxes [HJ, optical lattices Jl(| and electrons float- 
ing on liquid helium At present two-qubit gates 
were experimentally realized with cold atoms and 
the Grover algorithm has been performed for three qubits 
made from nuclear spins in a molecule [l9| . 

Thus there are two main lines in the present day quan- 
tum computer research: construction and development 
of efficient quantum algorithms and search for a optimal 
physical system with a future experimental realization of 
few coupled qubits. The first line has a strong mathc- 
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matical shade: indeed, it assumes that all qubits are per- 
fectly identical and the couplings between them can be 
operated also in a perfect way. The second line, in a large 
respect, is a part of experimental physics with few qubits. 
As a result, there is a broad open field for physical stud- 
ies of a realistic quantum computer with many qubits. 
Indeed, in reality the qubits are never perfect, the level 
spacing fluctuates from one qubit to another due to dif- 
ferent environment, there is also a residual interaction be- 
tween qubits which cannot be eliminated completely, the 
two-qubit gate operations are also not perfect. In prac- 
tice, at least around n = 1000 of such qubits are required 
to make a quantum computer more efficient than the ex- 
isting our days computers In addition to the above 
internal imperfections there is also decoherence due to 
coupling to external world which produces noise and dis- 
sipation. The effects of decoherence and two-qubit gates 
pulse broadening were numerically tested on Shor's algo- 
rithm 1 20 2l[] and were shown to play an important role. 
At the same time the estimates show that it is possible to 
have physical qubits with very long relaxation time, dur- 
ing which many gate operations can be realized p2[ . One 
of the most promising system looks to be nuclear spins 
in two-dimensional semiconductor structures [jll 12,23], 
However, the absence of external decoherence does 
not yet mean that the computer will operate properly. 
Indeed, internal imperfections with inter-qubit residual 
couplings J can strongly modify the ideal quantum reg- 
ister represented by noninteracting many-body (multi- 
qubit) states of ideal qubits. A simple estimate ]24| ] 
shows that the residual interaction J will be unavoid- 
ably much larger than the energy level spacing A„ be- 
tween adjacent eigenstates of a realistic/generic quan- 
tum computer . Let us assume that J is relatively weak 
comparing to one-qubit level spacing Ao. Then all Nu 
eigenenergies will be located in the energy band of size 
AE ~ uAq and the average multi-qubit level spacing is 
A„ m AE/N H ~ A 2~" < A . For the experimental 
proposals HJl2] we have Ao ~ 1 K so that for n = 1000, 
when Shor's algorithm becomes useful, the multi-qubit 
spacing is incredibly small A„ ~ 10 3 x 2~ 10 ' A ~ 10~ 298 
K. This value will definitely be much smaller than any 
physical residual interaction J. For the proposal 
with a distance between nuclear spins of r = 200 
and an effective Bohr radius of as = 30 A ( Eq.(2) 
of |l2|| ), the coupling between qubits (spin-spin interac- 
tion) is J ~ Ao ~ 1 K. By changing the electrostatic 
gate potential, an effective barrier between nuclei can 
be modified that can be viewed as a change of effec- 
tive electron mass possibly up to a factor of two. Since 
Joe (r/as) 5 ' 2 exp(— 2r/as)/as, and as is inversely pro- 
portional to the effective mass, this gives a minimal resid- 
ual spin-spin interaction of J ~ 10 -5 K ^S> A n . On the 
first glance this would lead to a natural/naive conclu- 
sion that at such residual interaction the ideal quantum 
computer eigenstates are strongly mixed and completely 



modified resulting in destruction of quantum computer 
hardware. In spite of this expectation it has been shown 
recently that the ideal qubit structure is much more ro- 
bust and in reality the quantum hardware melting and 
quantum chaos induced by inter-qubit interaction takes 
place at J > J c ~ Aq/u being exponentially larger than 
A„ p3]. This result for quantum chaos border in quan- 
tum computing is recently confirmed by more extended 
studies ]25| ] and opens a broad regime of parameters for 
which realization of a quantum computer is possible. For 
example, at Ao ~ 1 K and n = 1000 the critical coupling 
J c ~ 1 mK is compatible with the experimental proposal 
& 

The above result is closely related to the long term 
research of quantum many-body systems, started by 
Wigner |^(| interested in "the properties of the wave 
functions of quantum mechanical systems which are as- 
sumed to be so complicated that statistical considera- 
tion can be applied to them". As a result, the ran- 
dom matrix theory (RMT) has been developed to explain 
the generic properties of complex eigenstates and energy 
spectra of many-body interacting systems such as heavy 
nuclei, many electron atoms and molecules |27|. The suc- 
cess of RMT was so impressive |2f| that the conditions 
of its applicability to many-body systems were not really 
realized until recently. Indeed, for example, in nuclei the 
density of states grows exponentially with excitation en- 
ergy and at a first glance the RMT is valid as soon as 
the interaction matrix elements U become comparable to 
the level spacing between many-body states A„. How- 
ever, in nature we have only two-body interaction and 
therefore, while the size of the Hilbert space grows ex- 
ponentially with the number of particles n, the number 
of nonzero interaction induced matrix elements grows not 
faster than n 2 . To study the spectral statistics in this sit- 
uation a two-body random interaction model (TBRIM) 
was introduced and it was shown [^9|j3C| l that in the limit 
of strong interaction the RMT remains valid even if the 
full Hamiltonian matrix is exponentially sparsed. How- 
ever, much more time was needed to understand the case 
of relatively weak interaction and to find the critical U c 
above which quantum chaos and RMT set in. Contrary 
to naive expectation U ~ A„ assumed by many au- 
thors until recently (see e.g. |3l|-|33]|) it now became clear 
that the quantum chaos border is exponentially larger 
than A n since only transitions between directly coupled 
states, and not the bulk density of states, are important 
for level mixing. As a result the interaction should be 
larger than the energy level spacing A c between directly 
coupled states U > U c ~ A c ^> A„ to change the level 
spacing statistics P(s) from the Poisson distribution to 
the Wigner-Dyson one and to generate quantum ergod- 
icity and chaos for eigenstates. As for my knowledge, for 
the first time this condition for quantum chaos in many- 
body systems had been formulated by Aberg who also by 
numerical simulations for a nucleus model had found the 
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change in the level statistics across this border [Q . Due 
to that I will call this condition Aberg criterion. 

In spite of other Aberg's papers (e.g. his result 

was not broadly known for the community probably be- 
cause his studies were mainly addressed to nuclear physi- 
cists and the numerical results were not so confirma- 
tory at that time. Also in nuclei the interaction is usu- 
ally quite strong and the quantum chaos sets in rather 
quickly. In fact it was shown by numerical simulations 
of other models that in an isolated many-body Fermi 
system a sufficiently strong interaction can induce dy- 
namical thermalization with the Fermi-Dirac distribu- 
tion |3^j3^]. Independently , the line of research related 
to the problem of two interacting particles in a random 
potential J3(| showed that a two-body interaction can 
lead to a number of unexpected results and that disor- 
der/chaos can strongly enhance the interaction. The last 
property in fact had been known from the spectacular 
Sushkov-Flambaum enhancement of weak interaction in 
nuclei |37j. The analytical studies of few particles with 
two-body random interaction U |3§| l showed that the mix- 
ing of levels, quantum chaos and RMT statistics appear 
only if U > A2, where A2 > A3 > A4... are the level 
spacing between 2-,3-,4-particle levels respectively. This 
result was generalized for many particles and the quan- 
tum chaos border U c ~ A c was proposed and confirmed 
by extensive numerical studies of the TBRIM model jlO| , 
independently of Aberg's papers. As it will be seen from 
the present paper the knowledge obtained for quantum 
many-body systems can be successfully used for such a 
new direction of research as quantum computing. 

In this paper I review the recent developments of quan- 
tum chaos theory in many-body systems obtained in 
Toulouse and show their links and importance for a quan- 
tum computer which also represents a many-body system 
with exponentially large Hilbert space. The paper is or- 
ganized as follows. In the next Section the analytical and 
numerical results are presented for emergence of quantum 
chaos in many-body systems. In Section III the standard 
generic quantum computer (SGQC) model is described 
and the results of Section II are generalized for quantum 
chaos border for quantum computing. In Section IV the 
time evolution in the regime of quantum chaos is studied 
and different time scales imposed by chaos and decoher- 
ence on quantum computing are discussed. The paper 
ends by the concluding remarks in the last Section. 

II. ABERG CRITERION FOR EMERGENCE OF 
QUANTUM CHAOS IN MANY-BODY SYSTEMS 



particle states is A 2 ~ V/m 2 , A 3 ~ V/m 3 respectively 
with A ^> A2 3> A3. The two-body matrix elements 
written in the noninteracting eigenbasis are supposed to 
be random with a typical value U12 ~ U and U23 ~ U for 
interaction between lst/2d and 2d/3d particles respec- 
tively. At the same time the interaction matrix element 
Z7i3 between lst/3d is taken to be zero for simplicity (the 
final result remains the same for t/13 ~ U). For two in- 
teracting particles, e.g. lst/2d, the levels become mixed 
by interaction and RMT spectral statistics sets in when 
the interaction becomes larger than the level spacing be- 
tween two-particles states: U12 > A2. On the contrary 
for U12 <C A2 the perturbation theory is valid and the 
eigenstates with interaction are determined by one nonin- 
teracting eigenstate. In the case of three particles the sit- 
uation is more complicated since the three-particle states 
are not directly coupled by two-body interaction. 

The matrix element U3 between 3-particle levels can 
be found by the second order perturbation theory. In 
this way the matrix element between initial state 1 123 > 
and final state 1 1'2'3' > is given by diagram presented in 
Fig.l with intermediate state |1'23 >. 

3 ^{E 1 +E 2 + E 3 - E v - E- 2 - E 3 ) ~ A U 

It is important that the summation is carried out only 
over single particle states 2 and the sum is mainly deter- 
mined by a term with a minimal detuning in denominator 
being of order A. As a result the level mixing sets in for 
U3 > A3 that gives the quantum chaos border (3f|: 

U ~ U c ~ A 2 > A 3 (2) 

This means that the 3-particle levels are mixed only when 
the interaction mixes two-particle levels that is the conse- 
quence of the two-body nature of interaction. In a similar 
way for few particles n ~ 3 the border U c ~ A2. This 
conclusion was confirmed in [p9l. 




Following |38[ let us consider first a case of three parti- 
cles located on m one-particle orbitals with energy level 
spacing A ~ V/m. For simplicity we assume the parti- 
cles to be distinguishable that however is not of principle 
importance for m ^ 1. The level spacing between 2,3- 



FIG. 1. Diagram for the effective 3-particle matrix element 
U 3 in (@), after [38]. 

Let us now consider a more general case of TBRIM 
p9 pTJ with n Fermi particles distributed over m en- 
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ergy orbitals e m i, m' — l,2,...,m. These energies are 
randomly and homogeneously distributed in the inter- 
val [0, mA] with spacing A. The total number of many- 
body states is N = m!/(n!(m — n)\) and they are cou- 
pled by random two-body matrix elements which value 
is uniformly distributed in the interval [—U,U]. Due to 
the two-body nature of interaction the number of multi- 
particle states coupled by interaction, or the number of 
direct transitions, is K = 1 + n(m — n) + n(n — l)(m — 
n)(m — ri — l)/4 pi] ]. All these transitions occur inside a 
two-body energy interval B = (2m — 4) A around the en- 
ergy of initial multi-particle state. For large m and n the 
number of transitions K is much smaller than the size of 
the matrix N but is much larger than the number of dif- 
ferent two-body matrix elements N 2 « m 2 /2. The Fermi 
energy of the systems is ep ss nA and the level spacing 
in the middle of the total energy band is exponentially 
small A„ s» (m — n)nA/N. 



exponentially larger than the level spacing A„ between 
multi-particle states. 
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FIG. 2. Dependence of the rescaled critical interaction 
strength U c /B, above which P(s) becomes close to the 
Wigner-Dyson statistics, on the number of directly coupled 
states K for 4 < m < 80 and 1/40 < n/m < 1/2. The line 
shows the theory (|) with C = 0.58, after [40]. 

At the middle of spectrum the density of directly cou- 
pled states is p c = K/B. According to the usual pertur- 
bation theory, these states, if they would be alone, are 
mixed if the coupling U is larger than their level spacing 
A c = 1/pc- This is the Aberg criterion |j34] , ^5| , which 
was also independently proposed in f40(| . This relation 
is also in agreement with the arguments given above for 
few particles. As a result, the onset of quantum chaos is 
expected for U > U c : 



U r = CA r = C 



B 
K 



2C 
p 2 n 2 



(3) 



where C is a numerical constant to be found and p 2 ~ 
N2/B w m/4A is the two-particle density at m 3> n ^> 1. 
It is important to stress that the critical coupling U c is 




3.0 4.0 

Log( p p ) 



6.0 



FIG. 3. Dependence of the rescaled IPR £/U 2 on p c p n ' 
layer model data for n = 3, A = 1 and 40 < m < 130 (o); 
n — 4 and 30 < m < 60 (diamonds). The straight line gives 
theory (|5|). Insert shows U c /A vs. p c A in log- log scale for 
the same parameters; the straight line is the fit U c = 0.62/p c . 
After [43]. 

The most direct way to detect the emergence of quan- 
tum chaos is by the change in the probability distribu- 
tion P(s) of nearest level spacings s, where s is mea- 
sured in units of average spacing. Indeed, for inte- 
grable systems the levels are noncorrelated and char- 
acterized by the Poisson distribution Pp(s) — exp(— s), 
while in the quantum chaos regime the statistics is close 
to the Wigncr surmise Pw( s ) = (7rs/2) exp(— 7rs 2 /4) 
pd[ . To identify a transition from one limiting case 
to another it is convenient to introduce the parameter 
V = f*°(P(s)-Pw(s))ds/ f° a {Pp( S )-P w { S ))ds, where 
So = 0.4729... is the intersection point of Pp(s) and 
Pw(s)- In this way 77 varies from 1 (P(s) — Pp(s)) to 
(P(s) = Pw(s) ) and the critical value of U c can be de- 
termined by the condition r](U c ) = r\ c = 0.3 jib]] . In fact 
the choice of rj c affects only the numerical constant C in 
(^). The chosen r\ c — 0.3 is close to the value r\ c w 0.2 
for the critical statistics at the Anderson transition on 
a 3-dimensional disordered lattice |42]| . The results of 
extensive numerical studies, performed in [^0| , are pre- 
sented in Fig. 2 and confirm the Aberg criterion ([|) in a 
large parameter range with C = 0.58 |^3| . Of course, the 
data |fi"o|| are obtained for a much larger parameter range 
than in ]34|,[35|]. Nevertheless, the above C value is in 
a good agreement with a numerical factor found in |35| ] 
(see Eq. 22 there), which gives C « 0.7, if to take into 
account that C is defined via the same average values of 
square matrix element. A similar value is also found in 
more advanced studies for n — 3,4 in the layer model 
approximation with the states selected in the energy in- 
terval A (see insert in Fig. 3) [H . Also the above studies 
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pi] , [35 40 4~4| s how that contrary to the sharp Anderson 
transition 42 a smooth crossover from one statistics to 
another takes place at U ~ U c (see, however, and 
references therein). 

The Aberg criterion can be applied not only to ex- 
cited states but also to low energy excitations near the 
Fermi level. Suppose we have Fermi gas with a tem- 
perature T ep. Then, according to the Fermi-Dirac 
distribution, the number of effectively interacting parti- 
cles is Sn ~ Tn/eF with the density of these two-particle 
states p 2 ~ T/A 2 and p c ~ p 2 Sn 2 ~ A(T/A) 3 . The 
total excitation energy is SE ~ TSn ~ T 2 /A. As a 
result the interaction induced/dynamical thermalization 
and the quantum chaos set in |40| only for 

SE > SE ch » A(A/f/) 2 / 3 ; T > T ch w A(A/[/) 1 / 3 (4) 

These relations follow also from the Aberg papers [ ^4|j35| 
even if they were not written directly there. The numer- 
ical constant in corresponds to r\ c — 0.3 p5[ . Below 
the border (Q) the eigenstates are not ergodic and the 
interaction is too weak to thermalize the fermions even 
if the multi-particle level spacing is exponentially small 
(A„ oc exp(-2.5(J£;/A) 1 / 2 (47)). After §§§|Q the 
dependence (Q) was also obtained in p6| . 

In the quantum chaos regime U > U c the local den- 
sity of states is described by the Breit-Wigner distribu- 
tion with the energy width T given by the Fermi golden 
rule T = 2nU 2 p c /?> @]. The value of T determines the 
spreading width of eigenstates mixed by interaction. The 
number of noninteracting eigenstates contributing to a 
given eigenstate can be measured through the inverse 
participation ratio (IPR) £ = 1/Sil a j| 4 i where aj are 
probability amplitudes in the noninteracting eigenbasis. 
The mixing of all levels in the interval V gives §3 



spacings between the two states of one qubit. The Hamil- 
tonian of this model reads: 



£ w Tp n w 2U p c p„ 



(5) 



where the numerical factor is taken in analogy with the 
known result for band random matrices. This analytical 
relation is in a good agreement with the numerical data 
shown in Fig. 3. It is important to note that at U = U c 
exponentially many states are mixed by interaction. The 
width r has also another important physical meaning: it 
determines the chaotic time scale r x « 1/r after which 
an initial noninteracting eigenstate disintegrates over ex- 
ponentially many (£) eigenstates of interacting system 
(here and below h = 1). 



III. STANDARD GENERIC QUANTUM 
COMPUTER MODEL 



In p4j the standard generic quantum computer 
(SGQC) model was introduced to describe a system of n 
qubits containing imperfections which generate a resid- 
ual inter-qubit coupling and fluctuations in the energy 



(6) 



i<3 



where the a% are the Pauli matrices for the qubit i and 
the second sum runs over nearest-neighbor qubit pairs on 
a two-dimensional lattice with periodic boundary condi- 
tions applied. The energy spacing between the two states 
of a qubit is represented by Ti randomly and uniformly 
distributed in the interval [A - S/2,A + 6/2]. The 
detuning parameter 5 gives the width of the distribu- 
tion near the average value Ao and may vary from 
to Ao. Fluctuations in the values of Ti appear gener- 
ally as a result of imperfections, e.g. local magnetic field 
and density fluctuations in the experimental proposals 
Iprlp^l . The couplings Jjj represent the residual static 
interaction between qubits which is always present for 
reasons explained in the introduction. They can origi- 
nate from spin-exciton exchange [ll|Jl^| , Coulomb inter- 
action ||, dipole-dipole interaction j[L7| , etc... To catch 
the general features of the different proposals, Jy are cho- 
sen randomly and uniformly distributed in the interval 
[—J, J]. This SGQC model describes the quantum com- 
puter hardware, while the gate operation in time should 
include additional time-dependent terms in the Hamilto- 
nian (g) and will be studied separately. At J = the 
noninteracting eigenstates of the SGQC model can be 
presented as \ipi >= \a±, a n > where au = 0, 1 marks 
the polarization of each individual qubit. These are the 
ideal eigenstates of a quantum computer called quantum 
register states. For J 7^ 0, these states are no longer 
eigenstates of the Hamiltonian, and the new eigenstates 
are now linear combinations of different quantum register 
states. The term multi-qubit states is used to denote the 
eigenstates of the SGQC model with interaction but also 
for the case J = 0. It is important to stress that the quan- 
tum computer operates in the middle of energy spectrum 
where the density of states is exponentially large and it is 
natural that the quantum chaos initially sets in this bulk 
part of the spectrum. In this respect low energy exci- 
tations are not important contrary to fermionic systems 
discussed above. 

It is interesting to note that when one site in (^) is 
coupled with all other sites and S — 2Ao then the sys- 
tem becomes equivalent to the quantum version of the 
classical Shcrrington-Kirpatrick spin glass model stud- 
ied in |48| . For such a quantum spin glass shard in the 
middle of the spectrum the coupling matrix element is 
U ~ J, the spacing between directly coupled states is 
A c ~ Ao/n 2 , since each state is coupled to n(n — l)/2 
states in the energy band of the order of Ao. There- 
fore, according to ([|) the quantum chaos and ergodicity 
emerge for J > A /n 2 >• A„ as it was shown analytically 
and numerically in p8[. 
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A similar analysis can be done for the SGQC model 
(0). Indeed, for 6 <C A and J < 8 the total spectrum 
is composed of n bands with inter-band spacing 2Aq and 
band width -Jn8. Within one band one quantum register 
state is coupled to about n states in an energy window 
of 28 so that A c ~ S/n and the quantum chaos border is 
given by |3 



J c = CqS/n 



(J) 



where C q is a numerical factor 



All above arguments 
remain valid up to 5 — Ao when the bands become over- 
lapped. It is important that J c 3> A„ and that J c de- 
creases with S. The last property is natural since at small 
S the states in one band are quasi-degenerate and it is 
easier to mix levels. 




FIG. 4. Dependence of r\ on the rescaled coupling strength 
J I J c for the states in the middle of the energy band for 
n — 6(*),9(o), 1 2(triangles), I5(squares); 8 — Ao. The upper 
insert shows log(J c /Ao) (diamonds) and log(J cs /Ao) (trian- 
gles) versus log(n); the variation of the scaled multi-qubit 
spacing A n /Ao with log(n) is shown for comparison (+). 
Dashed line gives the theoretical formula (Q) with C q = 3.16; 
the solid line is J cs — 0.41Ao/n. The lower insert shows 
log(J cs /Ao) versus log(5/Ao) for n = 6(*),9 (o), 12 (trian- 
gles); straight lines have slope 1. After [24]. 

The direct numerical simulations for the quantum 



chaos border in quantum computer are done in 24 for 
the SGQC model. The change in the level spacing statis- 
tics P(s) in the band center with the growth of residual 
interaction J can be quantitatively characterized by the 
parameter ?y. To suppress fluctuations P(s) is obtained 
by averaging over 5 < Njj < 4 x 10 4 random realiza- 
tions of Ti,Jij so that the total spacing statistics was 
10 4 < Ns < 1.6 x 10 5 . Also P(s) is determined inside 
one of the symmetry classes of (|^) with odd or even num- 
ber of qubits up. The dependence of r] on J is presented 
in Fig. 4. The variation of critical coupling with n and 8 
can be determined from the condition ?y( J c ) = 0.3. The 
data obtained are in a good agreement with (Q) with 



C q = 3.16 and clearly show that J c ^> A„ (see inserts in 
Fig. 4). According to Fig. 4 the transition is sharp in 
the limit of large n in contrast to the smooth crossover 
in the TBRIM. This difference is due to local interac- 
tion between particles in (^) contrary to the long range 
interaction in the TBRIM. 

In the limit 8 <C A and J <C A the coupling 
between different energy bands is negligibly small. In 
this case to a good approximation the SGQC Hamilto- 
nian (|^) can be reduced to the renormalized Hamilto- 
nian Hp = S^+jPfci/Pfe where Pj~ is the projector on 
the k th band, so that qubits are coupled only inside one 
band. For an even n this band is centered exactly at 
E = 0, while for odd n there are two bands centered at 
E = ±A , and we will use the one at E = — A . Such a 
band corresponds to the highest density of states, and in 
a sense represents the quantum computer core. The de- 
pendence of critical coupling, determined via rj (J c ) = 0.3, 
on the number of qubits n is shown in Fig. 5 being in 
good agreement with the theoretical quantum chaos bor- 
der (0). It is important to note that at 8 = the pa- 
rameter 77 = and the eigenstates are chaotic |2jJ in 
agreement with (|t]). 
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FIG. 5. Dependence of log(J c /5) (diamonds) and 
log(J cs /<5) (triangles) versus log(n); the variation of the scaled 
multi-qubit spacing (log(A n /<5)) with log(n) is shown for com- 
parison (+). Dashed line gives the theoretical formula ([?]) 
with C q = 3.3; the solid line is J cs — 0.415/n; the dotted 
curve is drawn to guide the eye for (+). After [25]. 

The drastic change in the level spacing statistics is in 
fact related to a qualitative change in the quantum com- 
puted eigenstate structure. For J« J c the eigenstates 
are very close to noninteracting multi-qubit states \ipi >, 
while for J > J c each eigenstate \4> m > becomes a mix- 
ture of exponentially many \ipi >. It is convenient to 
characterize the complexity of an eigenstate \4> m > by the 



quantum eigenstate entropy S q 



J2 l W irn log 2 Wi, 



where Wi m is the quantum probability to find the state 
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\ipi > in the eigenstate \<j> m > (W im = \ < 4>i\(j> m > | 2 ). 
In this way S q = if \<j) m > is represented by one \ipi >, 
S q = 1 if two > with equal probability contribute to 
one \<p m > and the maximal value S q = n corresponds 
to a mixture of all Nh states in one \<p m >. A mix- 
ture of two states \tpi > is already sufficient to mod- 
ify strongly the quantum register computations and it is 
natural to determine the critical coupling J cs by the con- 
dition S q (J cs ) = 1. The results for J cs dependence on 
n and S are shown in Figs. 4,5. They clearly show that 
J cs w 0.13J C w OAS/n > A„. 

In the quantum chaos regime at J > J c one eigenstate 
is composed of £ states \ipi > mixed in the Breit-Wigner 
width T - J 2 /A c - J 2 n/S. As for the TBRIM the IPR 
is exponentially large and is given by (0) p9| . For J > S 
the interaction becomes too strong and T ~ Jy 7 "- ©• 
The pictorial image of the quantum computer melting is 
shown in Fig. 6. For J > J c the eigenstates become 
very complex and the quantum computer hardware and 
its operability are destroyed by residual inter-qubit inter- 
action. 




FIG. 6. The quantum computer melting induced by the 
coupling between qubits. Color represents the level of quan- 
tum eigenstate entropy S q , from bright red (S q ~ 11) to blue 
(S q — 0). Horizontal axis is the energy of the computer eigen- 
states counted from the ground state to the maximal energy 
(~ 2nAo). Vertical axis is the value of J/Ao, varying from 
to 0.5. Here n = 12, 8 — Ao, J c /Ao = 0.273, and one random 
realization of (0) is chosen. After [24]. 



IV. TIME SCALES FOR QUANTUM CHAOS 
AND DECOHERENCE IN QUANTUM 
COMPUTING 

The results of above section definitely show that the 
optimal regime for quantum computing corresponds to 



the integrable region below chaos border J < J c . How- 
ever, it is possible that for certain experimental propos- 
als it will be difficult to avoid the quantum chaos regime 
J > J c . Therefore, it is also important to understand 
during what time scale the quantum chaos becomes com- 
pletely developed. As for TBRIM this chaotic time scale 
t x is given by the decay rate T from one quantum register 
state to all others [§|||§: 

T x «i/r, r~j 2 n/s, r-j^ (8) 

where the expressions for T are given for J < 6 and J > 5 
respectively. For d — Aq a similar estimate for t x was 
given in |5Q| . 

The detailed numerical studies of chaotic disintegra- 
tion of an initial state \x(t = 0) >= \ipi >, correspond- 
ing to the quantum register state io, are done in [p5| . 
There for J > J c the behavior of the projection proba- 
bility F iio (t) — | < t()i\x(t) > | 2 is found and it is shown 
that the probability to stay at the initial state Fj j (t) 
decays rapidly to zero with the time scale r x w 1/r 
(see Fig. 14 there). The growth of the quantum en- 
tropy S(t) = - Y,i F iio (t) log 2 F Uo (t) with time is shown 
in Fig. 7. It clearly shows that after a finite time compa- 
rable with t x exponentially many states are excited and 
the computer operability is quickly destroyed. 
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FIG. 7. Time-dependence of the quantum entropy S(t) for 
J/ S — 0.4 > Jc/5 and n = 16 (diamonds), n = 15 (squares), 
n — 12 (triangles), n = 9 (circles), n = 6 (*). Average is made 
over 200 initial states io randomly chosen in the central band. 
Insert shows the same curves normalized to their maximal 
value. After [25]. 

Above we discussed the emergence of quantum chaos 
induced by inter-qubit coupling in an isolated quantum 
computer. It is possible to assume that this process can 
model to a certain extend the effects of decoherence in- 
duced by coupling to external world. Indeed, the two- 
body nature of interaction is also valid for external cou- 
pling. In addition it is also important to discuss another 



7 



type of danger for quantum computing which is not nec- 
essary related to mixing and complex structure of eigen- 
states. Indeed, the algorithms constructed for quantum 
computing, e.g. Shor's or Grover's algorithms, are based 
on ideal qubits which have 6 — 0. However, for 6 > 0, 
and even if J = 0, there is a phase difference A</> be- 
tween different states in one band which grows with time 
as Acf> ~ AE t, where AE is the energy difference be- 
tween states in the band and its maximal value can be 
estimated as AE ~ 6y/n. It is natural to assume that as 
soon as the accumulated phase difference becomes com- 
parable to 1 the computational errors become too large 
for correct quantum simulations. Therefore, in addition 
to the chaotic time scale t x , there exists the phase coher- 
ence time scale which is determined by the condition 
A0(t0) ~ AEt^ ~ 1. This scale is finite even at J = 
and can be estimated as ~ l/(6y/n) since the band 
width is AE ~ 6^/n and the computer usually operates 
with all states inside the band. For small J this estimate 
is still valid and in this case < t x for J < 6/n 1 ' 4 '. 
On the contrary for 6 = and J > the energy band 
width is AE ~ J\fn and according to (ih both scales are 
comparable ~ t x . While both time scales t x and 
are important for the quantum computer operability it is 
clear that the effects of quantum chaos on the scale t x are 
much more dangerous since after this scale an exponen- 
tially large number of quantum register states are mixed 
as it is shown in Fig. 7. On the contrary on the scale 
only the phases are changed but not the number of states. 
Therefore, it is natural to expect that phase spreading 
can be more easily suppressed by error-correcting codes 
than the onset of quantum chaos. 
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FIG. 8. Time-dependence of the energy E =< n 2 /2 > in 
the Chirikov standard map ([h]) in the chaotic regime with 
K = 5, k — 20 for 1000 orbits homogeneously distributed 
at n — for t = 0. The straight line shows the theoretical 
diffusion E = k 2 t/A. All velocities are inverted after t=150 
with computer accuracy 10~ 12 that destroys time reversibility 
after 20 iterations. After [53]. 



These results show that there are two important time 
scales t x and r^. The effect of external noise and deco- 
herence can be also characterized by two different scales. 
Indeed, from one side the noise gives some effective rate 
Tt with which a multi-qubit state decays to other states 
that determines the time scale r x ~ 1/Tt- It is clear that 
Tt oc n since the noise acts on all qubits independently. 
At the same time the noise gives some effective diffusion 
in energy which can be estimated as De = oj 2 Vt where 
lu is a typical energy change induced by noise during time 
1/Tt- For example, noisy fluctuations of Jy in (|^) give 
transitions with an energy change to ~ 6. As a result, 
AE ~ yU^> and Acj) - D e 1,2 t^/ 2 - 1. This deter- 
mines two time scales for external decoherence: 

r x «l/r T , r « 1/(lo 2 T t )^ 3 (9) 

where t x is related to the decay and relaxation, while 
determines the phase coherence. Let us note that the 
phase decoherence was discussed for electrons in metals 
at low temperature |3lf and was observed experimen- 
tally, see e.g. |p^| . However, there Tt was independent 
of number of electrons, while for quantum computing Tt 
is proportional to the number of qubits since the global 
coherence of multi-qubit states should be preserved. 




FIG. 9. Same as in Fig. 8 but for the corresponding quan- 
tum system ( p"l| ) with K = 5, k = 20, % = 1 realized on 
the same computer. The vertical line marks the inversion of 
time made by replacement tp — > tp* . The quantum dynamics 
remain perfectly reversible. After [53]. 

The time scales (^) play an important role for quan- 
tum computing. Indeed, the gate operations should be 
fast enough comparing to these scales to allow to realize 
error-correcting codes and to avoid a destruction of 
operability after the time scales (|^). Here, it is important 
to stress an important property of quantum evolution for 
which there is no exponential growth of errors, contrary 
to the classical dynamics. An illustration of this fact is 
based on the dynamics of the Chirikov standard map p3| : 
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n = n + ksm(6 + Tn/2) 1 d = d + T(n + n)/2 (10) 

Here, (n, 0) is a pair of momentum and phase vari- 
ables and bar denotes the new values of variables af- 
ter one period of perturbation. The classical dynam- 
ics depends only on the chaos parameter K = kT and 
for K > 0.9716.. the dynamics becomes globally diffu- 
sive in n. In this regime the classical trajectories are 
exponentially unstable and have positive Kolomogorov 
entropy. Due to that the computer errors grow exponen- 
tially quickly in time that in practice destroys the time- 
reversibility of the map ([l0]). This fact is illustrated on 
Fig. 8 |S^] where the iterations are done on a computer 
with errors of the order of 10~ 12 . Due to exponential in- 
stability of classical orbits the time-reversibility for orbits 
with inverted momentum (n — > — n) completely disap- 
pears after 20 iterations. The situation is absolutely dif- 
ferent for the corresponding quantum dynamics described 
by the unitary evolution operator for wave function tp on 
one iteration 54 1: 



tp = exp(-iTn 2 /4) exp{-ikcos9) exp(-iTn 2 /4) (11) 

where n = —id/d6 and K = 1. Here the quantum dy- 
namics simulated on the same computer remains exactly 
reversible as it is shown in Fig. 9. 

The above results clearly show that there is no expo- 
nential instability in quantum mechanics H] . This gives 
a direct indication that the errors in quantum computing 
can be efficiently corrected since e.g. Shor's algorithm re- 
quires relatively small number of gate operations Nq for 
number factorization: ~ 300L 3 where L is a number 
of digits in the number to be factorized |5). In a sense 
the Shor algorithm is exponentially fast while the errors 
grow only as a power of time that opens broad perspec- 
tives for quantum computing. This however requires a 
development of efficient error-correcting codes. There is 
also another general question related to the uncertainty 
relation between energy and time. Indeed, the quantum 
computing during time At allows to resolve energy levels 
only on the energy scale AE ~ 1 J At. This means that 
the resolution of exponentially small spacing AE ~ A„ 
requires exponentially long time analogous to the Heisen- 
berg time scale in the quantum chaos tn ~ 1/A n . Ap- 
parently only after this time scale all exponentially large 
information hidden in the Hilbert space can become avail- 
able. However, it is possible that on shorter time scales a 
useful information can be extracted in a way unaccessi- 
ble to classical computers, e.g. Shor's factorization. But 
a question about how many such exponentially efficient 
algorithms can be found remains open. 



V. CONCLUSION 

In this paper the conditions for emergence of quan- 
tum chaos, ergodicity and dynamical thermalization in 



many-body quantum systems with interaction and dis- 
order are presented. The Aberg criterion represents the 
main condition for onset of quantum chaos and different 
checks performed for various physical systems confirm its 
validity. The generalization of these results allows to de- 
termine the quantum chaos border for quantum comput- 
ing. In particular, it is shown that the critical coupling 
between qubits, which leads to quantum chaos and quan- 
tum computer hardware melting, drops only linearly with 
the number of qubits and is exponentially larger than 
the energy level spacing between eigenstates of quantum 
computer. In this sense the ideal multi-qubit structure 
is rather robust in respect to perturbations that opens 
broad possibilities for realization of quantum computers. 

Of course, the optimal regime for quantum computer 
operability corresponds to the integrable regime below 
the quantum chaos border. In this respect the quantum 
chaos is a negative effect for quantum computing which 
should be eliminated. Here it is possible to make an 
analogy with the development of classical chaos. This 
phenomenon also has negative effects for operability of 
plasma traps and accelerators. In fact the very first 
Chirikov resonance-overlap criterion had been invented 
in the pioneering work ]55[| for the explanation of exper- 
iments on plasma confinement in open magnetic traps 
and later found broad applications in accelerator physics 
p6[ . Since 1959 it is the only simple physical criterion 
which allows to determine the chaos border in classical 
nonlinear hamiltonian systems J55|,^6|,[53| and it merits to 
be marked by the Nobel symposium organizers. Indeed, 
the interest to the quantum chaos appeared only after a 
deep understanding of classical chaos had been reached. 
In a similar way it is possible to think that the deep 
understanding of quantum chaos in many-body systems 
in future will allow to quantum computers operate in a 
better way. 

The majority of results on quantum chaos in quantum 
computing presented in this paper were obtained with 
Bertrand Georgeot and it is my pleasure to thank him 
for the fruitful collaboration. I am also thankful to Oleg 
Sushkov for valuable discussions of various questions of 
interaction and disorder. I thank the IDRIS in Orsay and 
the CICT in Toulouse for access to their supercomput- 
ers. This research is partially done in the frame of EC 
program RTN1-1999-00400. 
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